Towards predictive modelling of near-edge structures in electron energy loss spectra 

of A1N based ternary alloys 
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Although electron energy loss near edge structure analysis provides a tool for experimentally 
probing unoccupied density of states, a detailed comparison with simulations is necessary in order 
to understand the origin of individual peaks. This paper presents a density functional theory 
based technique for predicting the N K-edge for ternary (quasi-binary) nitrogen alloys by adopting 
a core hole approach, a methodology that has been successful for binary nitride compounds. It is 
demonstrated that using the spectra of binary compounds for optimising the core hole charge (0.35 e 
for cubic Tii-^Al^N and 0.45 e for wurtzite Al^Gai-^N), the predicted spectra evolutions of the 
ternary alloys agree well with the experiments. The spectral features are subsequently discussed in 
terms of the electronic structure and bonding of the alloys. 

PACS numbers: 61.66.Dk, 71.15.Mb, 71. 20. Be, 71.20.Nr, 79.20.Uv, 81.15.Cd, 81.15.Gh 

Keywords: Density functional theory (DFT); Electron energy loss near edge structure (ELNES); Alloys; 
TiAIN; AlGaN 



I. INTRODUCTION 

Wurtzite aluminium nitride (w-AIN), gallium nitride 
(w-GaN) and their ternary alloy w-ALrGai-xN are im- 
portant materials for devices such as light emitting diodes 
(LEDs) and laser diodes (LDs)i. Cubic titanium ni- 
tride (c-TiN) and in particular its alloy with A1N, cu- 
bic Tii_ x Al x N, are widely used hard coatings due to 
their high hardness, corrosion and oxidation resistance 2 -. 
In both ternary alloys, a crucial requirement for getting 
the optimal application-tailored properties is an accurate 
control of both the composition and structure of the al- 
loy. 

Electron energy loss spectroscopy (EELS) is a pow- 
erful technique to microanalyse compositional features 3 -. 
The spectra can be recorded with very high spatial res- 
olution, thus taking advantage of the high sensitivity 
of EELS to local changes in the electronic structure of 
materials^ 5 -. Its subset, electron energy loss near edge 
structure (ELNES) reflects the density of unoccupied 
states and provides thus an experimental probe for this 
part of the electronic structure of materials. However, the 
interaction of high energy electrons with lattice atoms 
does not always have a straightforward interpretation. 
In order to understand the experimental data, measured 
EELS spectra need to be compared in detail with their 
calculated counterparts. 

X-ray absorbtion near-edge structure (XANES) is a 
closely related technique to ELNES. The edge shapes 
are very similar—, it also allows for studying polarisa- 
tion effects^, and it has been successfully applied to di- 



luted alloys^. Therefore, the here proposed methodology 
is highly relevant also for XANES. 

Recent developments of theoretical methods for solid 
state physics have provided the EELS community with 
increasingly reliable and comprehensive tools to simu- 
late ELNES. A particularly important example is the 
Telnes program which is distributed as a part of Wicn2k 9 , 
an all-electron full-potential linearised augmented plane- 
wave (FP-LAPW) code. It has been suggested in the 
literature^ that a calculation including a core hole pro- 
vides a better description of the excitation process by 
means of the standard (ground state) density functional 
theory (DFT). To create a core hole, one takes an elec- 
tron (or a fraction of it) from its ground state position (N 
ls-state in the case of N K-edge) and puts it in the lowest 
unoccupied state above the Fermi level (N 2p) or adds it 
as a background charge to keep the total charge of the cell 
neutral^. This can be easily realised in Wien2k where 
all electrons are accessible and explicitly treated. Re- 
cently, the core hole calculations have become routinely 
available also for pseudopotential code o n i 12 . 

A considerable effort has been spent on studying 
ELNES of binary nitrides with respect to their crystal 
structures (wurtzite, rock salt, zinc blende )^~— , polari- 
sation cffect a 13 ' 14 i 18 ~ — , doping 2 ^, and stoichiometry 2 - 2 -"— . 
The literature is vast and provides a good background for 
understanding the origin of peaks in ELNES in terms of 
bonding, and thus establishes a solid basis for a finger- 
print identification of materials and their properties. 

Nonetheless, there are only a few reports on the compo- 
sitional dependence of ELNES of ternary alloys. Keast 
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(a) B1 (NaCI structure) (b) B4 (ZnS structure) 




FIG. 1. Crystal structures investigated in this work: (a) cubic 
Bl (NaCI prototype), (b) wurtzite B4 (ZnS prototype), and 
(c) five-coordinated hexagonal Bk structure (BN prototype). 

et measured N K-cdge ELNES of In^Gai^N al- 

loys and correlated their findings with DFT calculations 
using, however, extremely small supercells. MacKenzie 
et al— used ordered structures and averaging of the 
boundary binaries spectra to get a guess on the shape of 
Tii^Al^N N K-edge. Holcc et al— showed that using 
small ordered cells failed to reproduce the experimental 
spectrum of an Al x Gai_ x N alloy. To date, a systematic 
study showing the effect of alloying on the ELNES shape 
( "evolution" of the edge) , as well as discussing the com- 
putational methodology in a recipe-like form, is missing. 
The present paper is aiming to fill this gap by facilitat- 
ing a semi- empirical approach in which the calculation 
parameters (e.g. the core hole charge) are first adjusted 
to reproduce the spectra of the boundary binary systems, 
and subsequently used to predict the N K-cdge evolution 
of ternary alloys. 

II. METHODOLOGY 
A. Calculation details 

The individual structures are modelled with super- 
cells constructed using a special quasi-random structure 
(SQS) approach 2 ^ All alloys considered in this paper 
are quasi-binary which means that mixing of elements 
(cither Ti and Al or Al and Ga) takes place only on one 
sublattice (bigger atoms in Fig. [TJ; the other sublatticc 
is fully occupied with N atoms. 3x3x2 (36 atoms) 
and 2x2x2 (32 atoms) supercells were used for the cu- 
bic Bl and wurtzite B4 modifications, respectively. The 
short range order parameters were optimised for pairs 
up to the fourth order, triplets up to the third order and 
quadruplets up to the second orde r 28 ' 29 . More details 
about the cells and the process of their generation can 
be found in Ref. H 

The cubic and wurtzite structures of the Tii-^Al^N 
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169 



TABLE I. Optimised lattice parameters and bulk moduli of 
the binary compounds used in this work. 

alloy were optimised with respect to volume (lattice con- 
stants) as well as to internal (local) relaxations, for which 
we used the VASP cod e 30 ' 31 together with the projector 
augmented wave pseudopotentials 3 ^. employing the gen- 
eralised gradient approximation (GGA) as parametrised 
by^ 3 -. We used 500 cV for the plane- wave cut-off energy 
and a minimum of ~ 600 /c-points-atom (usually more). 
Such parameters guarantee the calculation accuracy in 
the order of meV/atom. The obtained equilibrium lat- 
tice parameters (see Tab. [TJ) and energies are in good 
agreement with those previously published 3 ^— . 

TiN is mechanically unstable in the four co-ordinated 
wurtzite (B4) structure (Fig. [TJd) and relaxes into a five- 
coordinated Bk structure (Fig. Q}:). The reason for this 
is that the presence of d electrons favours a different hy- 
bridisation scheme (sp 3 d) than the tetrahedral sp 3 (see 
Refs. l27l and |37|) . Consequently, the wurtzite variant of 
the Tii-^Al^N alloy becomes unstable around x ps 0.6 
whereas for x J$ 0.6 the Bk structure is obtaine d 27 ' 38 . 
This is, however, the composition where also the phase 
transition to a lower energy, experimentally observed 
cubic variant happens 3 ^. Therefore, in the following 
we present only results for the high Al containing w- 
Tix^ALsN. 

As for the Al^Gai-^N alloy, we optimised the crys- 
tal lattices only for the boundary binary compounds 
(Tab. H|, A1N and GaN, and then used Vegard's rule to 
obtain the lattice constants for the intermediate compo- 
sition. This is justified by the work of Dridi et al— who 
showed for this alloy that the lattice parameters (unlike 
the band gap) exhibit a linear dependence on the com- 
position. 

Electronic properties and the ELNES spectra were cal- 
culated using the Wien2k code 9 - employing the GGA- 
PBE parametrisationiS, of the exchange-correlation po- 
tential. An equivalent of ps 900 fc-points within the 
whole first Brillouin zone of the unit cell, the expansion 
of the spherical harmonics up to I = 10 inside the non- 
overlapping muffin tin (MT) spheres, and -RMT^max = 7 
were usect^. The MT radii were automatically set by 
structGen (a part of the Wien2k package) to values 
« 1.70-1.80, ps 1.95-2.00, ps 1.85-1.95 and ps 1.90- 
1.95 a. u. for N, Ti, Al, and Ga atom, respectively. The 
spin polarisation effects were not taken into account. The 
core holes were implemented by reducing the N Is core 
level occupation on a specific site and putting the corre- 
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sponding charge in the background in order to keep the 
cell neutral. ELNES was calculated using the Telncs 
program, a part of Wien2k. This was repeated for all 
N sites in the supercell. The spectrum representing the 
particular alloy composition was obtained by averaging 
this set of N K-edges. 



B. Experiment 

In order to experimentally confirm the ab initio pre- 
dicted ELNES spectra, two material systems, Tii-^ALjN 
and ALjGai-zN were investigated. 

The Tii-zALjN samples were deposited in Leoben us- 
ing the plasma-assisted unbalanced magnetron sputter- 
ing technique^ 2 -. The variation of the Al mole fraction 
x in Tii„ x Al a: N was achieved by using powder metal- 
lurgically produced targets (PLANSEE AG, 99% pu- 
rity), with Ti/Al ratios of 1, 0.5 and 0.33, and man- 
ual placing of additional Ti or Al platelets (05 x 3 mm) 
in the racetrack of the targets, respectively. TEM sam- 
ple preparation was performed by Ar-ion thinning in a 
Gatan precision ion polishing system (PIPS) at 4 and 
2.2 kcV in plan view. The EELS measurements were car- 
ried out on a Cs corrected Jeol 2100F operated at 200 kV 
and equipped with a Gatan Tridicm GIF camera using 
nanobeam diffraction mode. This ensures high signal to 
noise ratios and makes it possible to acquire information 
from individual grains, necessary for the investigation 
of the polycrystalline Tii-^Al^N films. Thus, several 
different grains were measured for each alloy composi- 
tion to rule out possible orientation effects. The spectra 
were recorded with a dispersion of 0.3eV/channel and 
the energy resolution, measured by the full-width at half- 
maximum of the zero-loss peak, was 1.5-1.8 eV. The con- 
vergence and collection semi-angles during analysis were 
5 mrad and > 8 mrad, respectively. 

The Al^Gai-^N films were epitaxially grown at Cam- 
bridge by 6 x 2-inch Thomas Swan Close-Coupled Show- 
erhead metalorganic vapour-phase epitaxy (MOVPE) 
system^. A standard two-step growth method was used 
to deposit a 5 fim. thick GaN pseudo-substrate on (0001) 
sapphiroM on which the Al x Gai_ x N layers were grown 
at 1020 °C. The different compositions were obtained 
by varying the flow rate for Al and Ga precursors. The 
EELS spectra were obtained on a FEI Tecnai F20 mi- 
croscope equipped with a Schottky FEG source, Gatan 
Imaging Filter and operated at 200 kV. The electron 
beam was parallel to the (0001) direction. 

Prior to the comparison with the ab initio calculations, 
all measured spectra were corrected for the dark current 
and the channel-to-channcl gain variation. The pre-edge 
background was extrapolated using a power-law function 
and subtracted from the original data£. 



III. RESULTS 
A. Binary compounds 

The strategy adopted in this paper is to find calcula- 
tion parameters that reproduce the N K-edge ELNES for 
the binary compounds as closely as possible, and subse- 
quently use these settings for predicting the ELNES evo- 
lution of the alloys. Figure [2] shows how the edge shape of 
cubic and wurtzite/hexagonal A1N and TiN changes with 
increasing core hole charge from e (ground state) to 1 c 
(final excited state - full core hole) . The effect is stronger 
for A1N than for TiN. This is due to fast core hole screen- 
ing in TiN originating from its metallic character—. Nev- 
ertheless, some small changes can still be observed, e.g. 
the peak broadening for the cubic modification or the dis- 
appearing high-energy shoulder of the main peak for the 
hexagonal TiN with increasing core hole charge. Lazar 
et al— arrived at the same conclusion for c-TiN based 
on the comparison of their calculations with experimental 
measurements. 

The strong effect of the core hole charge on the A1N 
N K-edge shape has been discussed in the literatur o 12 ' 27 . 
A detailed analysis of the relative peak positions and in- 
tensities for the wurtzite A1N revealed that a core hole 
charge « 0.5-0.6 c reproduces the experimental ELNES 
the best^S. Such a comparison could still be misleading 
as the experimental results also depend strongly on the 
acquisition conditions 2 ^. Since the spectra around 0.5 e 
(Slater's transition state) for each allotrope look akin and 
are almost equally resembling the experiment, another 
criterion was adopted here. The edge onset, measured as 
the energy between the initial core state and the lowest 
unoccupied stated, is plotted as a function of the core hole 
charge in Fig. [3J This dependence is rather strong, and 
using the experimental value for the edge onset allows 
an optimal value of the core hole charge to be estimated. 
Taking 402 eV for w-AlNi^ (or this work) and 397 eV 
for c-TiNi&il (or this work) yields 0.45 e and 0.35 e, re- 
spectively, which are the values used in this work. 

B. Tii -xAl^N alloy 

The calculated evolution of the N K-edge for 
Tii-^ALjN is shown by solid lines in Fig. [4] for the cubic 
and wurtzite modifications. The raw ELNES was broad- 
ened with a Gaussian having 1 eV FWHM. Moreover, 
the curves of pure A1N [x = 1) were shifted along the en- 
ergy axis (as labelled in Fig. [4j to account for the abrupt 
change in the Fermi energy due to the development of 
the band gap (no Ti d-states present). 

The compositional step used for the cubic alloy is Ax = 
0.167 (Fig. HJi) ■ Three developments of the main peaks 
are predicted: (i) the double-maxima at 0-5 eV above 
Ep disappears with increasing A1N content, (ii) the peak 
at « 13 eV for TiN gradually moves to « 9eV for A1N 
and at the same time its intensity increases, and (iii) a 
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FIG. 2. Calculated N K-edge ELNES as a function of the core hole charge changing from a ground state calculation (Oe) to a 
full core hole (le): (a) cubic Bl A1N, (b) wurtzite B4 A1N, (c) cubic Bl TiN, and (d) hexagonal B k TiN. 
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FIG. 3. N K-edge onset (energy difference between the initial 
core state and lowest unoccupied state) as a function of the 
core hole charge for the two allotropes of each A1N and TiN. 
The inset shows zoomed-in region around the experimentally 
observed N K-edge onset energy. 



small hump at « 32 eV for TiN broadens to an almost 
undetectable background at x = 0.5. At the same time, 
a small hump develops at w 26 cV with increasing A1N 
content. The origins of these composition- induced peak 
variations can be tracked down to the changes in bonding 
in the alloy as is discussed later in section HVBI 

The situation is more complicated for the w- 
Ti^^Al^N alloy (Fig. [Dd. compositional step Ax — 
0.125). The calculated spectra suggest that the char- 



acteristic triple-peak character of the w-AIN N K-edge 
is levelled out with the addition of only 0.125 mole frac- 
tion of TiN. Additionally, the spectra do not show any 
clear trends in the peak development as in the cubic 
case, apart from the peak at « 20 eV, whose position 
is not influenced by the composition. This is most likely 
connected with local relaxations taking place similar to 
those reported for Nbi-^Al^N^. In particular, N sites 
near Al tend to have the four-coordinated (wurtzite-like) 
neighbourhood, while in the vicinity of Ti atoms, five- 
coordinated local neighbourhoods are preferred (hexago- 
nal Bk-like) . The structure is therefore much more sensi- 
tive to the actual arrangement of atoms in the supercell 
which is reflected, e.g. in much bigger scatter of the "op- 
timised" lattice constants for different arrangements of 
atoms in the SQS (with a constant composition 



To compare the calculated and measured N K-edge 
evolutions (Fig. a larger spectrometer broadening 
parameter of 1.5 eV was used. As a consequence, the 
double-maximum of the first peak in c-Tii-^Al^N at 
~ 0-5 eV above the Ep "smears out" and the measured 
shape is obtained. It is therefore concluded that the fine 
double-maximum character is not resolved due to ex- 
perimental limitations. The experimental spectra were 
smoothed and normalised to fit the intensity of the high- 
est peak of the simulated pattern in each individual case; 
no other operation was performed on them. The theo- 
retical spectra were, on the other hand, shifted by the 
calculated energy of the core-holed core level. The spec- 
tra thus obtained exhibit a very good correlation between 
experiment and theory. 
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Energy above the Fermi level [eV] Energy above the Fermi level [eV] Electron energy loss [eV] 

FIG. 4. Calculated N K-edge evolution in the (a) cubic and (b) wurtzite phase of Tii-^ALjN alloy. The spectra of pure A1N 
(x — 1) are shifted in order to correct for the opened band gap. The dashed line corresponds to a linear interpolation between 
the binary compounds spectra (see section llV A[) . (c) Comparison of the calculated (solid black line) and experimental N K-edge 
ELNES of Tii-^ALrN. Note that the structure changes from cubic (x < 0.7) to wurtzite (x > 0.7). 



C. Al^Gai-zN alloy IV. DISCUSSION 

A. Shape and evolution of the N K-edge ELNES 

The spectra of the binary TiN, A1N and GaN systems 
have been extensively discussed in the literature both 
from the experimental and theoretical perspective. De- 
spite that, several issues remain unclear, in particular the 
edge onset energy: its value varies in the literature within 
the range of several eV o 12 ' 16 ' 17 ' 21 ' 26 ' 27 ' 45 . Consequently, 
we used our own measurements to calibrate the calcula- 
tions instead of taking spectra from the literature. The 
lincshapes of binary w-AIN, c-TiN, and w-GaN resemble 
those previously published for the same materials. 

The edge evolution for c-Ti 1 _ 2 .Al ;1 ,N exhibits the same 
trends as the one previously published by MacKcnzic 
et ai— (which, however, provided only one intermediate 
composition) . Gago et ai— used XANES to measure the 
N K-edge of Tii^A^N experimentally. Their XANES 
spectra have all the main features of our experimental 
as well as calculated N K-edge ELNES. Also in the case 
of w-ALjGai-j-N, the calculated evolution of N K-edge 
ELNES correlates with the here presented experimental 
data as well with those published previousl y 2 ^ 5152 . 

MacKenzie et ai— tried to model the edge with a lin- 
ear interpolation of the boundary binary spectra. Having 
in mind the problems with the accurate edge onset deter- 
mination and the lack of cubic A1N (in the Bl structure) 



As another example, the semiconducting wurtzite solid 
solution of A1N and GaN is chosen to demonstrate the 
ability of the current approach to predict ELNES. In 
contrast to the meta-stable Tii-^Al^N, the wurtzite 
Al^Gai^N mixture is stable in this modification for all 
concentrations x. Therefore, no local distortions as in the 
case of the w-Tii-^Al^N alloy are expected, which results 
in a gradual N K-edge evolution as shown in Fig. [5^,. The 
intensities of the first and third peaks of the triple-peak 
shape characteristic for w-AIN decrease with decreasing 
A1N mole fraction, and become shoulders around a cen- 
tral peak, a shape typical for GaN. The peak at « 24 eV 
above the Ep in the w-AIN spectrum gradually moves to 
~ 26 eV for w-GaN with increasing GaN content. The 
edge onset moves towards the Fermi level reflecting the 
narrowing band gap from 4.2 eV (A1N) to 1.7 eV (GaN)^. 
This gradual tranformation of the spectra shape is traced 
down to the changes in electronic structure, see section 
ITVBl The predicted evolution of the N K-edge is again 
confirmed by the experimental observations (Fig. [SJd) . 
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FIG. 5. (a) Calculated evolution of N K-edge of the wurtzite ALjGai-zN alloy and (b) a comparison of calculated (solid black 
line) with experimental N K-edge ELNES shapes. The dashed lines correspond to linearly interpolated boundary spectra of 
A1N (x = 1) and GaN (x = 0). 



for getting a reliable binary spectrum, this approach is 
questionable. To demonstrate this further, Figs. [4^ and 
[S£i include the linear interpolations of the binary spectra 
(dashed lines). Although this may serve as a first (and 
quick) guess on what the evolution should look like, in 
many cases the relative intensities and/or positions of the 
peaks are not predicted correctly. 

Craven— showed using several binary transition metal 
nitrides (TMN) that the spacing between the double- 
maximum peaks increases with increasing lattice param- 
eter. This is not predicted for the c-Tii-^Al^N al- 
loy (Fig. |4|) where the lattice parameter decreases from 
4.25 A for TiN (x = 0) to 4.07 A for A1N (a; = 1)24, but 
the peak spacing is almost unaffected. The reason for 
this is that the bonding of various TMN is similar and the 
peaks follow small shifts of the density of states associ- 
ated with the varying valence configuration. On the con- 
trary, the peak shifts in the Tii-^Al^N evolution result 
in the first place from the changing character of bond- 
ing (see section IIVB[) . Consequently, when the ELNES 
of the two boundary binary systems are similar, the sim- 
ple approach of interpolating between the binary ELNES 
spectra^ is expected to give acceptable results, see e.g. 
Tii^V^N^ or Ini.^Ga^Ni^. The extremely small cells 
(1 In and 1 Ga atom for Ino.5Gao.5N) used in the latter 
reference satisfactorily reproduced the N K-edge evolu- 
tion, and a much more computationally expensive ap- 
proach (using the averaging of several N core-holed sites 
in SQSs as in the present paper) employed by Holec— is 



not necessary. 

On the other hand, when the spectrum evolution is 
pronounced, small ordered structures do not provide re- 
liable predictions. This has been shown by MacKcnzie 
et at— for the case of c-Tii-^ALjN and by Holec et at— 
for w-Al^Gai-^N. In such cases, the approach adopted 
here is essential. 



B. Electronic origin of the peaks 

There exists extensive literature on the origin of peaks 
for semiconducting III-N binaries. As summarised by 
Mizoguchi— using the overlap population analysis, the 
main peak structure (up to « 10 eV above the edge onset) 
reflects the anti-bonding N-cation interaction while the 
later peak (20-30 eV above the edge onset) corresponds 
to cation-cation (mostly anti-bonding) interactions. The 
difference between the A1N and GaN ELNES shapes can 
be traced down to the presence of the valence d-electrons 
in GaN which cause (slight) redistribution of the valence 
density of states, and consequently also the unoccupied 
density of states. The electronic structure of the valence 
band of InN is similar to that of GaN thus resulting in a 
similar ELNES spectrum. 

The ground state projected density of states (PDOS) 
in Fig. [5] helps to understand the meaning and evolution 
of peaks in ELNES of the c-Tii^Al^N alloy. The site 
and symmetry projected DOS were obtained by averag- 
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(a) x=0 (b) x=0.167 (c) x=0.333 (d) x=0.5 (e) x=0.667 




Energy above Fermi level [eV] 

FIG. 6. PDOS/ELNES spectra for o-Tii-^AlsN with (a) x = 0, (b) x = 0.167, (c) x = 0.333, (d) x = 0.5, and (e) x = 0.667. 
The peaks in N K-edge ELNES closely reflect the distribution of the unoccupied N p states. The overlap of the individual 
PDOSes suggest that the double peak at « 2-4 eV above Ef arises from hybridisation of N p and Ti d states while the peak 
at ~ 10 eV can be ascribed primarily to N p states hybridised with Al p states. 



ing the corresponding quantities over all sites occupied 
with the same specie. The final states of the N K-edge 
transition are unoccupied N p-states which clearly cor- 
relate with the ELNES. Based on the PDOS overlaps 
it can be concluded that the double-maximum structure 
just above the edge onset arises from the N p-Ti d-states 
interaction which agrees with the findings of Tsujimoto 
et al.— and Lazar et al.— for binary c-TiN. The second 
peak at around 10-12 eV above Ep have the strongest 
contribution from the N p-Al p-states interaction, the 
only exception being pure TiN (no Al present) where a 
small peak in Ti d-states at the same position can be 
detected. The different interactions contributing to this 
delayed peak are responsible for a sharper maximum with 
clear shoulders in the case of pure TiN while resulting in 
a rather broad (and nearly symmetric) shape when Al 
is present, see Fig. [6] The peak position changes by al- 
most leV upon adding x = 0.167 mole fraction of A1N 
to TiN, while further increase of A1N content results in 
only a small shift of the peak (0.5 eV for increasing x 
from 0.167 to 0.667). This can serve as an example why 
the simple interpolation between the properties of binary 
compounds (as suggested by MacKenzie et al.—) does 
not work. 

The bonding of the c-Tii-^ALjN alloy consists of a 
mixture of covalent and ionic type. The bonding of TiN 
has been discussed many times in the literature^. It 
was concluded that the covalent part is realised by the 
sp 3 d 2 orbitals (interaction of Ti 4s and 3d-states with N 
2p-states). Additionally, the interaction between Ti 3d 
orbitals with the t^ g symmetry causes a non-zero DOS at 
the Fermi level resulting in the metallic character of the 
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FIG. 7. Bader analysis of a charge transfer (absolute values) 
on atoms of individual species. 



compound. The influence of Al on bonding in the alloy 
was discussed by several author a 34 ' 35 ' 56 , generally show- 
ing a gradual weakening of the sp 3 d 2 hybridisation (which 
is reflected by the decreased intensity of the first double- 
maximum peak in the N K-edge ELNES, see Figs. Q] and 
[6]). Additionally, the Bader analysis^ as implemented 
in Wien2k shows that there is a significantly increased 
charge transfer from metallic sites to N resulting in a 
stronger ionic bonding with increasing AIN mole fraction 
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(see Fig. [7]). It is interesting to note that the different AI 
sites "provide" on average always almost the same charge 
to be transferred on N sites (practically no scatter around 
the mean values shown by the triangles in Fig.[7|), but the 
charge transferred from Ti sites is much more influenced 
by the alloy composition. This is likely to be due to dif- 
ferent degrees of hybridisation between Ti and N atoms 
depending on the actual neighbourhood of N atoms (i.e. 
second-order neighbours of Ti sites) . 



C. Influence of the local environment 

There is some controversy in the literature on how big 
the supercells should be in order to suppress the mu- 
tual interactions between core holes. For example, Mi- 
zoguchi et ai— and Tanaka and Mizoguchi— claimed 
that cells with more than 100 atoms are needed while 
32 atom cells were found sufficient by Lazar et ai— for 
GaN and by Holec et ai— for w-AlN. To address this is- 
sue we plotted the individual spectra resulting from core 
hole being placed on various N sites, and sorted them 
according to the number of the nearest neighbours of 
each specie (in total 4 for the tetrahedrally coordinated 
wurtzite structure (Fig. [8^) and 6 for the octahcdrally 
coordinated cubic structure (Fig. [8)3, c)) surrounding the 
particular N site with the core hole. The numbers of 
spectra corresponding to individual local environments 
(i.e. the nearest neighbour configuration) results from 
their real counts in the used supercells. Although con- 
figurations of the nearest neighbours of some N sites are 
the same, the higher-order nearest neighbour configura- 
tions differ which is why small variations between indi- 
vidual spectra labelled with the same local environment 
are obtained. The thick lines on top of each panels in 
Fig. [8] were obtained by averaging all the curves under- 
neath (their number is the same as the number of N sites 
in the supercell) , and thus account for the statistical dis- 
tribution of various local environments of N atoms. 

The graphs clearly demonstrate the huge differences 
between spectra depending on the local environment of 
the N site where the excitation takes place. At the 
same time one can see that almost doubling the num- 
ber of atoms in the supercell (from 36 in Fig. [8b to 64 in 
Fig. \8jp) does not alter the N K-edge significantly. The 
small changes might be due to having insufficiently big 
cell in the case of 36 atoms, but also could be due to 
a non-representative (i.e. not SQS-like) cell for the big- 
ger structure which is suggested, e.g. by not having the 
2A1, 4Ti local environment present. In summary, Fig. [5] 
demonstrates that (i) the local environment influences 
the final shape of the edge much more drastically than the 
actual cell size (provided the cell is big enough to model 
the "randomness" of an alloy), and (ii) the similarity of 
the curves from A1-, Ti-, Ga-rich local neighbourhoods 
to those of pure A1N, TiN and GaN, respectively, gives 
some grounds for the interpolation approach (see section 



D. Energy of the edge onset 

The edge onset energy is another important feature 
of the edge; for example, Mizoguchi et ai— predicted 
a range of ~ 2eV for the N K-edge onset of A1N de- 
pending on the crystal structure thus having a potential 
to distinguish between these allotropes. It is, however, 
not straightforward to define the edge energy, in partic- 
ular due to the ambiguity in the background subtraction 
as well as due to the background noise itself. Conse- 
quently, we have chosen the energy of the first inflection 
point above the edge threshold for the comparison be- 
tween experiment and theory (the edge onset is about 
2-4 eV below). 

It is not surprising that we obtained excellent agree- 
ment for the binary systems (< 0.1 eV for GaN and 
w 0.35 eV for A1N) since for these systems the edge onset 
energy was used as a fitting parameter for the core hole 
charge (see section UlI Aj) . However, both the experiment 
and the theory suggest that the energy of the inflection 
point does not vary too much with the composition. The 
variations within 0.2 eV (theory) and 0.4 eV (experiment) 
can be regarded as the accuracy of the present approach 
(due to, e.g. the background subtraction on the experi- 
mental side or the supercell design/size on the theoretical 
side) . 

A common method calculating the excitation energy 
is to calculate the difference between total energies of 
the initial ground state and the final state with a full 
core hol o 15 ' 46 ' 58 which, in principle, follows the excitation 
process^. This results in values of 384.5 eV and 368.2 eV 
for w-AIN and c-TiN, respectively, which are hugely un- 
derestimated as compared with the experimental values 
of 402 eV and 397 eV. The reason for this is that the ex- 
cited electron was put as the background charge (Fig. [SJd) 
rather than in the unoccupied states (Fig. [9J:). When 
excited to the unoccupied states values of 406.4 eV and 
404.6 eV for the w-AIN and c-TiN are obtained which 
are much closer to the experimental values. Rashkova 
et ai— showed that a further improvement (towards the 
experimental values) could be obtained by performing 
spin polarised calculations. Nevertheless, the edge shapes 
as well as the energies of the initial core levels are almost 
identical using both approaches (background charge vs. 
valence band) thus yielding comparable results (except 
for the edge onset), which is in agreement with Hebert— . 

When the edge onset is calculated as the total energy 
difference between the ground and full core hole states, 
its value is a given number without any degree of free- 
dom for adjustments. The corresponding ELNES shape 
then should be calculated with exactly 0.5 e core hole^i. 
This could be useful when estimating, e.g. ELNES of 
experimentally inaccessible phases. However, when one 
uses the core hole charge as a fitting parameter (as in 
this paper) then it is well justified that also the edge 
onset is not a unique number but instead a function of 
the core hole charge. For evaluation of the energy differ- 
ence between the core state and the lowest unoccupied 
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FIG. 8. The resulting spectra representing the alloys (thick curves on top of each panel) and spectra calculated by placing the 
core hole on individual N sites and sorted accorting to the nearest neighbour configurations, (a) Wurtzite Alo.5Gao.5N, (b) 
cubic Tio.5Alo.5N supercells with 36 atoms, and (c) cubic Tio.5Alo.5N supercells with 64 atoms. 
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Fermi 
level 




FIG. 9. (a) A ground state and an excited state with a full core hole: (b) the excited electron is put as the "background charge" 
(i.e. effectively removed from the system) while (c) it is put in the valence band (into prior the excitation unoccupied states). 
The scheme corresponds to a metallic system. 



state, however, the approach with background charge is 
more appropriate since, for example, in the case of the 
conductive TiN it allows to get the energy of the lowest 
(originally) unoccupied state (compare Figs. [9Jd and c). 



V. CONCLUSIONS 

In this paper we have demonstrated a semi-empirical 
approach for predictive calculations of the N K-edge 
ELNES of various classes of alloys (cubic vs. wurtzite, 
metallic vs. semiconducting). We used fractional core 
holes with charges, carefully adjusted according to the 
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edge shape and the onset energy, to reproduce experi- 
mental spectra. Subsequently we utilised these to model 
the ELNES spectra of alloys using the special quasi- 
random supercells. We introduced core holes on all 
individual N sites, and by averaging these spectra we 
achieved a representative alloy spectrum. A comparison 
with the experimental measurements (for c-Tii-^ALjN 
and w-Al x Gai_a;N systems) yielded an excellent agree- 
ment on both the edge shapes (including peak positions 
and relative intensities) as well as the edge onset ener- 
gies. Finally, we related the individual peaks in the N 
K-edge ELNES to various interactions between cations 
and N atoms demonstrating that the decrease in inten- 



sity of the N K-cdgc structure 3 eV above the edge on- 
set reflects a weakening of the sp 3 d 2 hybridisation with 
increasing Al content in Tii-xAl^N. 
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